Distribution function of mesoscopic hopping conductance 
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We study by computer simulation distribution functions (DF) of mesoscopic hopping conductance. 
The DFs obtained for one-dimensional systems were found to be quite close to the predictions of the 
theory by Raikh and Ruzin. For D — 2, the DFs both for narrow system and thin film look similar 
(and close to the ID case) .The distribution function for the conductance of the square sample is 
nearly Gaussian. 
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Mesoscopic conductance fluctuations in the insulating 
regime of small, disordered transistors were first observed 
by Pepper Q in GaAs MESFETs and then studied in de- 
tail in Si MOSFETs by Fowler, Webb and coworkers ^ in 
the early 1980s. Extremely strong random fluctuations, 
spanning several orders of magnitude, were observed at 
low temperatures in the conductances of narrow-channel 
devices as the gate voltage was varied. The explanation 
was provided by Lee |^ who proposed a model in which 
electrons move by variable-range hopping (VRH) along 
a one-dimensional (ID) chain. A number of elementary 
hopping resistances, each depending exponentially on the 
separation and energy difference between sites, are added 
in series to give the overall resistance of the chain. In this 
model it is assumed that, because of the extremely broad 
distribution of the elementary resistors, the total chain 
resistance can be well approximated by that of the single 
most resistive hop. The fluctuations then arise as a con- 
sequence of switching between the pairs of localized sites 
responsible for the critical hop as each elementary resis- 
tance reacts differently to a change in the chemical po- 
tential. These fluctuations are therefore of "geometrical" 
origin, arising from the random positioning of localized 
sites in energy and space, as distinct from the "quan- 
tum" nature of the tunneling mechanism which would 
be strongly affected for example by an applied magnetic 
field. Serota, Kalia and Lee Q went on to simulate the 
ensemble distribution of the total chain resistance R and 
its dependence on the temperature T and the sample 
length L. In their ensemble, the random impurities are 
distributed uniformly in energy and position along the 
chain. In experiments a single device is generally used, 
so that the impurity configuration is fixed, and fluctua- 
tions are observed as a function of some variable external 
parameter such as the chemical potential. An ergodicity 
hypothesis is then invoked to the effect that the same 
ensemble is sampled in both cases, something that has 
been verified experimentally by Orlov et al. [1). Using 



the natural logarithm of the resistance, the authors of 
Ref. [Q obtained for the mean and standard deviation: 
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where ^ is the localization radius and Tq is the character- 
istic temperature for Mott VRH: Tq — l/ksp^ {p is the 
density of states at the Fermi energy) . It can be seen that 
the size s of the fluctuations decreases extremely slowly 
with length, a result characteristic of ID which was flrst 
pointed out by Kurkijarvi The explanation is simply 
that exceptionally large resistance elements, even though 
they may be statistically rare, dominate the overall re- 
sistance since they cannot be by-passed in this geometry. 
The averaging assumed in the derivation of Mott's hop- 
ping law for ID does not occur and the total resistance 
takes on the activated form of the largest individual ele- 
ment. 

A detailed analytical treatment of this model was un- 
dertaken by Raikh and Ruzin (RR) who divided the 
problem into a number of length regimes. Their theory 
introduces the concept of the "optimal break" , the type 
of gap between localized states (on an energy versus po- 
sition plot) which is most likely to determine the overall 
resistance. The optimal shape of such a state-free region 
has maximal resistance for the smallest area and turns 
out to be a rhombus. A sufficiently long chain will have 
many such breaks in series to give a most probable resis- 
tance 
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where Rq is the prefactor in the Mott VRH formula, and 
p — R/Rq. This formula is valid only for substantially 
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long samples. More exactly, it is valid when i/ ^ 1, where 
z/ is a parameter defined implicitly by 
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When the expected number of optimal breaks in the 
chain becomes of order one {v < 1), which corresponds to 
the normal experimental situation, the resistance of the 
chain is determined by a few sub-optimal breaks of which 
the expected number occurring in a chain of length L is 
approximately one. In this case, the most probable value 
of the resistance (or its logarithm Q = In = In p) is 
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The probability distribution function (DF) for the 
quantity f{Q) is best written in terms of v and a new 
parameter A 



T ' 



(6) 



For < 1 it is given by the following integral: 
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where /(A) is a function with a peak close to A = and 
width determined by i^. There is a simple relationship 
between v and the variance of Q: 
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This theory is equally applicable |g[| to the case of the 
transverse conductance ct of a thin film or barrier. In- 
stead of a sum of series resistances, the required quantity 
is the sum of parallel conductances representing conduct- 
ing chains of hops traversing the film. Whereas in ID 
the total resistance is determined by the blocking effect 
of the critical hop, here the total conductance is domi- 
nated by an optimal "puncture" : an uncommonly high- 
conductance hopping chain through the barrier which ef- 
fectively shorts out all other current paths. On a loga- 
rithmic scale, since Inp = —In a, the DFs for the two 
geometries are simply reflections of each other. The vari- 
ation of the width and peak position with ^ and p is dif- 
ferent, however, in the two cases. In ID the importance 
of blocking resistors adds weight to the contribution of 
extremely high resistances and produces a long tail out to 
low values of Inc. For a short 2D barrier the DF has the 
opposite asymmetry with a tail out to high conductances, 
reflecting the effect of punctures in shorting out less con- 
ductive paths. In fact the form of the DF is universal, the 



theory requiring only that the elementary quantities to 
be summed are independent and come from an exponen- 
tially wide distribution. The microscopic details of the 
conduction mechanism enter only into the dependence of 
v and A on external parameters such as the temperature 
and magnetic field. The requirement of independence in 
the case of the barrier means that conductive chains must 
be sufficiently far apart, which should be satisfied for a 
barrier with a sufficiently large aspect ratio W/L. We use 
the description "short 2D" for this short-, wide-channel 
geometry to distinguish it from the square 2D geometry 
in which conduction is via an interconnected percolation 
network. 

The aim of the present work is to do numerical sim- 
ulation on the ID and 2D mesoscopic systems in Mott 
hopping regime, and find the distribution functions of 
conductance in these systems. We start by replacing 
the transport problem with a random-resistor network 
in which the hopping between sites i and j is equivalent 
to having a resistor pij such that 

\np,j = 2ad+{\E, - fi\ + \Ej - p\ + \E, - Ej\)/2kT, 

(9) 

Here pij is the resistor between sites i and j, a is the 
inverse localization length, d is the distance of two local- 
ized sites, Ei and Ej are energies of site and p, is the 
chemical potential, T is the temperature. Thus the meso- 
scopic system is reduced to a random resistor network. 
To find the resistance of the network, the resistors joining 
electrodes are selected in ascending order until the first 
percolation path connects the reservoirs. The resistance 
of the percolation path is taken to be the resistance of 
the entire system. 

In the simulation, first, a number of different impu- 
rity configuration is generated. Sites are randomly dis- 
tributed along dimensions of the system and their en- 
ergies are chosen randomly from a uniform distribution 
between —0.5 ~ -fO.5. For each configuration we ran- 
domly chose the position of the chemical potential p. 
Thus we can consider the chemical potential distributions 
(for a fixed impurity configuration) and the ensemble and 
chemical potential distributions. Typical results for a 
given configuration and given value of chemical potential 
are presented below. For a ID system oi L — 1000 When 
chemical potential p — and temperature T = 0.001, by 
the method described above, the threshold is In pc — 11.9. 
Among about 500,000 resistors, there are 53 resistors sat- 
isfying unequality Inpy < Inpc, among which there are 
12 resistors in the percolation path. The profile of their 
resistances is shown on Fig. la. We see that the sin- 
gle largest hop in fact determines the conductance of the 
system. With temperature increasing, the resistors in 
percolation path become closer in resistance, which can 
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been seen on Fig. lb, which presents the results of similar 
calculations, but with T = 0.01. 
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FIG. 1. The ratio of individual resistances in the percola- 
tion path to the path resistance: a)T = 0.001; b)r = 0.01. 
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FIG. 3. The conductance of ID system for v = 3.6: a) 
ensemble distribution function; b) chemical potential distri- 
bution function. The histograms in both figures are the nu- 
merical result, the solid lines correspond to the fitting curves, 
and the dash-dot lines are the prediction by the RR theory. 



We simulate the ID system at first, and compare the 
numerical results with the RR theory. The ID system 
has a length of 1000, and a localization length of 50. We 
consider three cases: 

1. T = 0.001, which gives u = 0.225; 

2. T = 0.01, which gives u = 3.6; 

3. T = 0.015, which gives u = 5.8. 

For each temperature we consider 1000 ensembles, and 
chemical potential range is = —0.1 ^ -1-0.1. The re- 
sults of the first case is shown below. 
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FIG. 2. The conductance of ID system for f = 0.225: a) 
ensemble distribution function; b) chemical potential distri- 
bution function. The histograms in both figures are the nu- 
merical result, the solid lines correspond to the fitting curves, 
and the dash-dot lines are the prediction by the RR theory. 



The numerical simulation results for = 3.6, where 
according to the theory, the distribution function should 
be Gaussian, are presented below. 



The result of T = 0.015 leads to the similar figures 
above except for thiner and higher DFs, which shifts a 
little towards the higher conductance. 



For 2D system, we consider three particiilar cases: nar- 
row 2D system, square sample and thin film. We expect 
that the DF of narrow 2D system should be close to that 
of ID system and according to the RR theory, the thin 
film DF on a logarithmic scale should be a reflection of 
the ID case distribution function. It is natural to expect 
that the DF of the square system would be close to a 
Gaussian. For a 2D system all parameters are chosen to 
be the same as in the ID case, except w = 100 in the case 
of narrow sample; w = 1000 and L = 100 for thin film 
and w — L — 1000 for square 2D system. The number 
of impurity configurations = 50 in each case. The results 
are shown below. 
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FIG. 4. The ensemble distribution function for the conduc- 
tance of 2D system: a) thin film; b) square sample; c) long 
system. Histograms are numerical results. The solid curve in 
b) is the Gaussian fitting, in c) is the RR theoretical fitting 
for the ID case. 

We see that the situations of narrow 2D is close to the 
ID case as expected, and the normal 2D is close to Gaus- 
sian. But, unexpectedly, the DF for thin film is similar 
to the ID case distribution function and not a mirror 
reflection of it, as argued by RR. 

In conclusion, the paper has studied the DF of the 
conductance in mesoscopic systems by numerical simula- 
tion. Wo have found that the distributions obtained by 
choosing randomly the chemical potentials (for a fixed 
impurity configuration), which corresponds to a typi- 
cal experimental situation, coincide with those obtained 
when both impurity configuration and chemical poten- 
tial is chosen randomly, in agreement with the ergodicity 
hypothesis. The DFs obtained for one-dimensional sys- 
tems were found to be quite close to the predictions of 



the theory by Raikh and Ruzin. For D = 2, the DF 
both for narrow system and thin film looks similar (and 
close to the ID case). The distribution function for the 
conductance of the square sample is close to Gaussian. 
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